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1 Introduction 

The role of control allocation (CA) in modern aerospace vehicles is to compute a command 
vector 6 C € m na that corresponding to commanded or desired body-frame torques (momcn s) 

r = [ L M IV ] T to the vehicle, compensating for and/or responding to inaccuracies in off-line 
nominal control allocation calculations, actuator failures and/or degradations (reduced cffcctiv- 
ncss), or actuator limitations (ratc/position saturation). The command vector 6 C may govern c 
behavior of, c.g., acrosurfaccs, reaction thrusters, engine gimbals and/or thrust vectoring. Typi- 
cally, the individual moments generated in response to each of the n a commands docs not lie stric y 
in the roll, pitch, or yaw axes, and so a common practice is to group or gang actuators so that a 
one-to-one mapping from torque commands r c to actuator commands 6 C may be achieved m an 

off-line computed CA function. n . . , . 

We shall assume the existence of an off-line computed nominal linear affine CA function 


= F{x)t c + do(z) 


( 1 - 1 ) 


where r c is the commanded torque vector, x is a vehicle state vector, d 0 » » tnm (neutral torque) 
vector and Fix) is a matrix of nominal control allocation gains. One may interpret the columns of 
Fix) as a set set of gains defining “ganged” actuators for each control axis. The resulting nominal 
autopilot/control allocation block diagram is shown in Figure 1. The vector n in Figure 1 refers 
to the body torques induced on the vehicle by the actuators. Ideally, the control allocation matrix 
F{x) would be chosen to be the pseudo-inverse G{x)' of the Jacobian matrix 

G(x) = [ TSjf ] | s € m 3 *"" 

where n„ is Ihc number of actuators: that is, we wish to design the control allocation matrix F(x) 
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Figure 1: Nominal autopilot/control allocation block diagram. 


such that 


G(i)F(i) = 


1 0 0 
0 1 o 
0 0 1 


( 1 . 2 ) 


so that the induced body-frame vehicle torques r 6 match the commanded body Hmvcva, 

due to sensor inaccuracy, modeling errors, and data compression in the alloc^ f ™ C 
ideal condition (1.2) cannot be achieved. As such, semidefimte programming techniques [HSZOlj, 
[BEFB94] arc used to design control allocation matrices F(x) that adueve 

||G(z)F(x) - I\\ < l-7(*) 

for some positive constant 7 to achieve the “best possible” allocation given off-line data^ 

The linear affine nominal control allocation law (1.1) by itself is inadequate for the control 

control allocation problem for four reasons: 

1. It fails to respond to torque allocation errors that can be detected on-line. 

2. It fails to take into account saturation issues. 

3. It fails to respond to on-line detected failures in actuators, and 

4. It fails to provide a framework to work with discrete-valued (on-off) actuators such as reaction 
thrusters. 

Dynamic control allocation [HC02] may be used to compensate for torque allocation <«o»dct«^ 
on-line. DCA treats torque allocation error as a unknown additive uncertainty in the system 

Jacobians . , . 

Gir ue{ x ) — C{x) + Acf(x) 

where Ac; is an unknown gain that respects the condition 

II Gjrue {x)F (ar) - /|| < 1 

over all operating conditions. If this condition is not met, 

olicatcd requiring the use of on-line system identification [CPM95] [CBP98], [CG86J, [HJN91J. 
£t3d™ tnators may be dealt with by either osing pul^width modnlafon to emnlate 
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continuous valued actuators or, where actuator structure will not permit the use of PWM, one may 
use these actuators as “back-up” to the continuous valued actuators as in [HSZ01]. We address in 
this paper items 2 and 3 above by autonomous rcconfigurablc control allocation (ARCA). 

An optimal solution of a constrained control allocation problem involves the solution of convex 
programming problems [Buf97], [Dur93], [Enn98]. Some heuristic approximate allocation solutions 
arc presented in [BD95] such as computation of the attainable moment subset, generalized inverses, 
and daisy chaining. Alternatively, adaptive control may be used in tandem with linear system 
theory techniques in an attempt to avoid and/or compensate for actuator saturation [CGD 98]. In 
this paper we build on the work of [BD95], [Dur93], [BP98], [BD95], and [PBOO] to present an on- 
line autonomous rcconfigurablc control allocation technique that is computationally tractable and 
practical for use in closed loop with a robust autopilot (attitude control law) such as sliding mode 
control [SMJ+98], [BLM99], [SHJOO]. Our technique makes use of a fast quadratic programming 
iteration step so that the actuator command 5 C results in a vehicle body torque n that tracks the 
commanded torques r c when they lie within the attainable moment set or else approximates the 
commanded torques in a least-squares sense when they arc not in the attainable moment set (the 
underlying quadratic programming problem is infeasible). 

2 Autonomous Reconfigurable Control Allocation (ARCA) 

Rcconfigurablc control refers to the ability of a control allocation law to continue to maintain 
tracking of the moment command r c in the face of actuator failurcs/dcgradation. In the short 
term, actuator saturation is indistinguishable from actuator failure, since in both eases an additional 
constraint is entered into the control allocation law. However, in the ease of an actuator failure, 
the constraint is permanent and thus requires coordinated treatment between the control allocation 
law, the attitude control law, and the guidance law. 


2.1 Problem statement 


We formally describe the ARCA problem as follows. We shall denote the sequence of autopilot 
torque commands as r c {k) and the sequence of CA generated actuator command vectors as 6 c (k). 
In order to accomodate actuator rate limits, we shall compute the actuator command vectors 
recursively, i.c., 

SAk) = 6 c {k - 1) + S e (k) 


where S e {k) is an update to the previous actuator command vector <5 c (fc - 1). We shall omit the 
dependence on the time index k where it is clear by context. We shall make use of the following 
definitions: 


Definition 2.1 The actuator status vector 5 slai {t) has entries in the range of [0, 1] where <W,» = 0 
reflects complete failure (actuator i has no impact on vehicle body torques) and £ siaij = 1 reflects 
nonminal operation of actuator i . 


Remark 2.1 We shall assume that the the actuator status vector S, lal is made available to the 
CA module by, c.g., a vehicle health monitoring system or an on-line system identification module. 
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Definition 2.2 The effective actuator Jacobian G(x) is the system Jacobian 

GW = [ & ] 

evaluated under current operating conditions, including effects of actuator degradation and failure. 

□ 

Remark 2.2 Notice that the effective actuator Jacobian G(x) is unknown prior to flight time. If 
the nominal system Jacobians G{x) « F{x)' arc available, then G{x) = G(x)diag(J» to t). 

Definition 2.3 Given a sampling interval T, the current actuator position 6 and vectors of actuator 
maximum position values 5+, minimum actuator position values 6 , and maximum actuator rates 
<Wx> the next- step feasible set Aj{5 c ) is defined as 

A f{6 c ) = |tf e : S~ < 5 C + 5 e < 6 + and |J e - <*l < T^max} 

Remark 2.3 The set A f {5{kT)) is the set of legal next-step actuator commands S c {{k + 1)T) given 
current actuator position values S[k). 

Denoting J c = <5 c (fc - 1), 5 e (k), the rcconfigurablc control allocation problem can thus be ex- 
pressed as the quadratic programming problem [Luc84] 

min (<5 C + 8 C ) T W\ (J c + ^e) + 8 e T W2^e 
subject to G(x)5 e = At 

and 5 e e Ay(5 c ) 

where Wi is a position penalty on the command vector 8 c {k), W 2 is a rate penalty on the change 
in command vector 5 e , and At is a torque command update, usually (but not always) defined as 
A T (fc) = T c (k) - T c (k - 1). Exceptions to this rule arc discussed below. 

In the ease where the quadratic programming problem (2.1) is feasible, then the mmimization 
searches for an update 5 e that matches the commanded torque while reducing the magnitude 
(cost) of the control command. Conversely, if the ideal torque command is not feasible, i.c., no 
solution 5 e € A f {5 c ) exists to the equality constraint G{x)5 e = At, then it is necessary to relax 
the equality constraint and instead find a vertex (extreme point [Luc84]) of the feasible set A/(A C ) 
that minimizes the norm \\T err (k)\\ where T e ^ G(x)8 e - At. In this ease, the next value of At is 
selected to reflect both the update to the torque command t c and the “unallocated torque T e from 

the previous iteration. . - . 

Brute force application of standard quadratic programming techniques may not be desirable m 

the ARCA problem for the following reasons. 

1. Computation time: solution of the quadratic programming problem can requires several it- 
erations, each requiring solution of a linear system of equations Ax = b of dimension up to 
n a x n a . 

2 Repetition: the quadratic programming problem to be solved in our application is solved 
repeatedly for problems and conditions that do not vary greatly from one sample time to the 
next. For example, during times of highly aggressive maneuvers, it is likely that the same 
actuators arc saturated from one time step to the next. The similarity of these problems is 
not exploited by standard quadratic programming techniques. 
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3 Singular dual solution: because the number of constraints is always 4 times the number of 
unknowns, the dual problem [LucG9] used to solve for the corresponding Lagrange multiplier 
in a non-negative least squares problem is singular, which requires at the least an increase m 
the dimension of the problem solution via the use of “slack” variables. 

2.2 Fast QP solution to ARCA 

Given the concerns raised above, we propose the use of a fast QP solution method that differs from 
standard primal problem QP solvers in two ways: 

1. The solution (active constraint set [Luc84]) of the previous problem is used as a starting point 
of the current problem. 

2. The number of QP iterations is limited so that computational burden is reduced. 

We consider here the equivalent QP problem using the notation of [Luc84]. 

min J{x) = x T Qx + c T x (2*2) 

X 

subject to Ax = b 

x~ < x < x+ 

The unknowns x in problem (2.2) correspond to the actuator command deviations 6 e in the ARCA 
algorithm description. We shall proceed on the following assumptions: 

1. x~ < 0 < x + : i.c., the previous actuator command vector satisfies actuator constraint limits. 

This assumption may be violated in the ease of, c.g., initial transients due to engine failure, 
where additional constraints on differential thrust can be imposed upon failure. Nevertheless 
the above assumption will be reasonable in all eases where additional constraints arc not 
suddenly imposed on actuator command behavior. 

2. The torque command issued by the autocommandcr is feasible. This condition is can be met 
through the use of on-line computation of a local attainable moment set [•]. 

Based on the above assumptions, we may solve the QP problem (2.2) as illustrated in Figure 
2 The initial point x(0) by assumption satisfies the inequality constraints x < x(0) * , u 

may not satisfy the equality constraint Ax = b due to variations in the (effective) Jacobian and 
torque command from one time step to the next. The point w{0) is computed that minimizes the 

optimization 

min tz;(0) (2*3) 

w(0) 

subject to x4u;(0) = b 

with A = A, b=b. As illustrated in the figure, u;(0) will satisfy the equality constraints, but 
may not satisfy the inequality constraints. We therefore choose our next solution value x(l) - (1 
a)x(0) +ma(0) where w{0) is chosen so that x(l) lies at the edge of the feasible set for the inequality 
constraints. We then append a row to A and b corresponding to the new “active constraint, 
x, = x7 in the ease of the figure, and solve the minimization (2.3) again, now constrained to the 
intersection of the hypcrplancs Ax = b and *1 = x^ This process continues until the iterate w(k) 
is in the feasible set of problem (2.2). (Such a solution exists by assumption.) 
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Figure 2: Quadratic programming problem (2.2) solution method for the ARCA problem. The 
feasible solution set is denoted by the rectangle on the interior of the hypcrplanc Ax - b. 

Remark 2.4 Because the system Jacobians arc continuous (except in the ease of sudden failure), 
the active constraints from one iteration to the next will likely not diangc often As a r ^> 
above algorithm may be easily modified to solve the minimization (2.3) at most once per contr 
step while evaluating revisions to the active constraint set by projecting the gradient 


onto the null space of A 


^- = 2Qx + c 
ax 


■ = (j - A* A) (2Qx + c) 


and comparing the signs of the resulting gradient search direction * with the active constraints 
imposed on x from the previous iteration. 

3 Conclusions 

Our initial tests of our fast QP solution method on artificially generated problems arejery en- 
couraging. For the final version of this paper we shall present closed-loop results using the ARCA 
algorithm in closed loop with a high-fidelity model of the X-33 experimental launch vehicle. 
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